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c3^ '■ ABSTRACT 
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This paper presents a photometric and spectroscopic study of the bright blue 

' r-|'' echpsing binary LMC-SCl-105, selected from the OGLE catalog as a candidate 

host of very massive stars (> 30 Mq). The system is found to be a double- 



Dh' 



P \ lined spectroscopic binary, which indeed contains massive stars. The masses 



\-i 



Vh ; and radii of the components are Mi = 30.9 ± 1.0 Mq, M2 = 13.0 ± 0.7 Mq, 

and Ri = 15.1 ± 0.2 Rq, R2 = 11.9 ± 0.2 Rq, respectively. The less massive 



(N 



X 



star is found to be filling its Roche lobe, indicating the system has undergone 



K>- \ mass-transfer. The spectra of LMC-SCl-105 display the Struve-Sahade effect, 

T^ij- \ with the He I lines of the secondary appearing stronger when it is receding and 

f^ \ causing the spectral types to change with phase (08-1-08 to 07-f08.5). This 

|.^ \ effect could be related to the mass-transfer in this system. To date, accurate (< 

O \ 10%) fundamental parameters have only been measured for 15 stars with masses 

Q \ greater than 30 Mq, with the reported measurements contributing valuable data 
on the fundamental parameters of very massive stars at low metallicity. The 



results of this work demonstrate that the strategy of targeting the brightest blue 



H ' stars in echpsing binaries is an effective way of studying very massive stars. 

Subject headings: binaries: eclipsing - binaries: spectroscopic - stars: funda- 
mental parameters, individual (OGLE053448. 26-694236. 4) - galaxies: individual 
(LMC) 



^ Based on observations obtained with the 2.5 meter DuPont and 6.5 meter Magellan Clay Telescope 
located at Las Campanas Observatory, Chile. 

^Vera Rubin Fellow 
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1. Introduction 

The fundamental parameters of very massive stars (> 30 Mq) remain uncertain, de- 
spite the large impac t massive star s have in astrophysics, both individually and collec- 
tively (see review by iMasseyl |2003| ). The equations of stellar structure allow for stars 
with arbitrarily l arge masses, however th e mechanisms to form massive stars (accretion 



and r nergers; e.g. iBallv fc Zinneckerl |2005| ) and the associated instabilities (see lElmegreen 



2000l : IZinnecker fc Yorkd 120071 . and references therein) are not well understood, hindering 
theoretical predictions on the existence of an upper limit on the stellar mass. The "mass 
discrepancy" problem, i.e. the disagreement between masses derived from parameters de- 
ter nnn^d__by_^ttiiigstellarat^^ to spectra and from evolutionary tracks (see 



e.g. iRepolust et al.l 12004 : iMassey et al.ll2005l . for a comparison), still affects studies of single 
massive stars, even though s ignificant progress has been made i n both stellar atmospher e 
(see review by iHerrerd 120071 ) and stellar evolution models (e.g. iMeynet fc Maederl |2003| ). 
The parameters of single stars also suffer from suspected multiplicity, which in many cases 
cannot be deter mined. Pismis 24-1 demonstrates this problem: its inferred evolutionary 
mass > 200 Mq ( IWalborn et al.ll2002l ) contradicted the up per stellar rnass limit of ~ 150 M f;^ 
suggested by statistical arg uments based on observations (lFigerll2005l : lOey &: Clarkd l2005l ) . 
Maiz Apellaniz et al.l (120071 ) resolved it into a visual binary with the Hubble Space Telescope, 
thereby removing the discrepancy. One of its components is also a spectroscopic binary, il- 
lustrating the systematic effects often accompanying "single" stars. 

The only model-independent way to obtain accurate fundamental parameters of distant 
massive st ars and to resol ve the "mass discrepancy" problem is to use eclipsing binaries (see 
review by I AndersenI 1 1 99 ll ) . In particular, double-lined spectroscopic binary systems exhibit- 
ing eclipses in their light curves are extremely powerful tools for measuring masses and radii 
of stars. Specifically, the light curve provides the orbital period, inclination, eccentricity, the 
fractional radii and flux ratio of the two stars. The radial velocity semi-amplitudes deter- 
mine the mass ratio; the individual masses can be solved for by using Kepler's third law. 
Furthermore, by fitting synthetic spectra to the observed ones, one can infer th e effective tem- 



perat ures of the stars, solve for their luminosities and derive the distance (e.g. iBonanos et al. 



20061 ). The most massive stars measured in eclipsing binaries are galactic Wolf -Rayet stars of 



iOi 



Schnurr et al. 



WN6 ha spectral type: NGC3603-A1 (Mi = 116 ± 31 Mq, M2 = 89 ± 16 M^ 

2008h. and WR 20a (Mi = 83.0 ± 5.0 Mq and M2 = 82.0 ± 5.0 Mq IrIuw et al. 



2004 



Bonanos et al.ll2004l ) in Wester l und 2, presenting a challenge for both ste . 



massive star formation models (JYungelson et al.ll2008l : IZinnecker fc Yorke 



lar evolution and 



20071) and raising 



the issue of the frequency and origi n of "binary twins" (jPinsonneault fc Stanekll2006l : iLucy 
20061 : iKrumholz &: Thompson! 120071 ). Such systems are of particular inter est, since massive 



binaries might be progenitors of gamma-ra v bursts fe.g. iFrver etal.1 120071 ). especially in the 



case of Population III, metal-free stars (see lBromm fc Loebll2006l ). 
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Analogs of these heavyweight champions, if not more massive binaries, are bound to 
exist in the young massive clusters at the center of the Galaxy (Center, Arches, Quintuplet), 
in nearby super star clusters (e.g. Westerlundl, R136), in Local Group galaxies (e.g. LMC, 
SMC, M31, M33) and beyond (e.g. M81, M83, NGC 2403). A systematic wide-ranging 
survey of these clusters and galaxies is currently underway. The goal is to provide data with 
which to test star formation theories, stellar atmosphere and stellar evolution models for 
both single and binary stars as a function of metallicity, and the theoretical predictions on 
the upper limit of the stellar mass. The adopted strategy involves two steps: a variability 
survey to discover eclipsing binaries in these massive clusters and nearby galaxies, which is 
followed by spectroscopy to derive parameters of the brightest - thus most luminous and 
massive - blue systems. However, characterizing massive stars requires the availability of 
8-m class telescopes and high resolution near-infrared spectrographs (since massive stars in 
the Galaxy are ext incted an d extra galactic ones are faint) and has only become feasible in 
the past few years. iBonanosI (120071 ) demonstrated that this method efficiently finds massive 
candidates, by performing the first variability survey of the Westerlund 1 super star cluster 
and discovering 4 massive eclipsing binary systems. 

Figure [T] illustrates the extent of our knowledge of precise fundamental parameters of 
massive stars. It presents published mass-radius measurements from eclipsing binaries, ac- 
curate to better than 10 % for the more mass ive component. The zero-age main sequences 
(ZAMS) at both Z=0.02 JSchaller et al.lll992h and Z=0.008 JSchaerer et al.lll993h are over- 



plotte d as a reference. The Galactic data are mainly taken from the compilat i ons o 



Gonzalez et a 



(Il99l[) and ICiesI (l2003f) with additions and updates from IVitrichenko et al.l (12007), bu t also 



J2005h for the Large M agellanic Cloud (L MC), Iflarries et al.l J2003h and 



Andersen 



Hilditch et all (l2005f ) for the SMC, and iRibas et all (l2005f ) for M31. A literature search 
was done to include all accurate measurements of stars in eclipsing binaries with masses 
> 30 Mq, which are presented in Table [H This Table is, to my knowledge, complete at 
present and consists of only 14 very massive stars with better than 10% mass-radius m ea- 
surements, located in 3 galaxies. Of these, WR 20a and M33 X-7 ( lOrosz et al.l 120071 ) are 
the most massive and noteworthy. M33 X-7 contains a very massive 70.0 ± 6.9 M© 0-type 
giant and a record-breaking 15.65 M© black hole, challenging current evolutionary models, 
which fail to explain such a large black hole mass. Without accurate measurements for a 
large sample of massive stars, theoretical models will remain unconstrained. 

A survey to determine accurate parameters for several massive eclipsing binaries in the 
low metallicity (Z= 0.008) LMC was undertaken, with the purpose of increasing the sample 
and improving our understanding of these rare systems. Several candidates were s elected 
from the OGLE-II catalog of eclipsing binaries in the LMC (jWyrzykowski et al.ll2003l ) as the 
brightest systems with 5- 1/ < 0. LMC-SCl-105, or OGLE053448.26-694236.4, is one of the 
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brightest eclipsing binaries with Imax = 13.04 mag, Vmax = 12.97 mag, Bmax = 12.81 mag 
and a preliminary semi-detached classification. This work presents the analysis of the follow- 
up observations obtained for LMC-SCl-105. The paper is organized as follows: §2 describes 
the spectroscopy and data reduction, §3 the spectral classification, §4 the radial velocity 
curve, §5 the light curve analysis, §6 the evolutionary status, and §7 the conclusion. 



2. Spectroscopy 

A total of 9 spectra of LMC-SCl-105 near quadrature phases were acquired over 4 
runs on 2 telescopes at Las Campanas Observatory, Chile. In December 2005, January and 
February 2006, spectra were obtained with the Echelle spectrograph on the 2.5-m DuPont 
telescope. The 1" x 4" slit resulted in a spectral resolution of 10 km s~^ {R = 30000) or 
0.15 A at 4500A, as measured from the full width half maximum of the comparison lamp 
lines. In January 2006 and October 200 7, additional spectra were obtained with the blue 
and red sides of the MIKE spectrograph (JBernstein et al.ll2003l ) on the 6.5-m Magellan Clay 
telescope. The 1.0" x 5.0" slit used in 2006 resulted in a spectral resolution of 9 km s~^ 
{R = 32000) or 0.14A at 4500A, as measured from the full width half maximum of the 
comparison lamp lines. The 0.7" x 5.0" slit used in 2007 resulted in a spectral resolution of 
7 km s~^ {R = 41000) or O.llA at 4500A. Table [2] summarizes the log of the observations, 
specifying the telescope and instrument used, the exposure times and final signal to noise 
(S/N) ratio per pixel measured on the merged spectra. 

Th e Magellan spectra were extracte d using the MIKE reduction pipeline written by D. 
Kelson (JKelson et al.ll2000l : lKelsonll2003l ). The extracted orders for each star were averaged, 
normalized and merged. The wavelength coverage of the final merged spectra is 3900 — 5050A 
in the blue and 5040 — 7150A in the red. The DuPont spectra were reduced, extracted and 
wavelength calibrated using the noao.imred. echelle package in IRAFn. Cosmic rays were 
removed from the two dimensional images with the algorithm of iPychI ( 120041 ). The orders 
were finally averaged, normalized and merged, yielding a wavelength range 3700 — 9000A. 
The heliocentric radial velocity corrections for each star were computed with the IRAF 
rvsao.bcvcorr routine and taken into account in the subsequent radial velocity determination. 
Narrow nebular emission, likely originating from a nearby HII region, is present in the 
spectra and was removed. Finally, the Ca II H and K and Na I D lines exhibit multiple 
absorption features, corresponding to Galactic and LMC interstellar material with radial 
velocities between and 315 km s~^. 



^IRAF is distributed by the NOAO, which are operated by the Association of Universities for Research 
in Astronomy, Inc., mider cooperative agreement with the NSF. 
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3. Struve-Sahade Effect and Spectral Classification 



Inspection of the quadratu re spectra reveals that LMC-SCl-105 exhibits the "Struve- 
Sahade effect" (ISticklandlll997l ). This term describes the variable strength of the spectral 
lines of the secondary star (or primary star in some cases) in a double-lined spectroscopic 
binary (see lHowarth et al.lll997l . and references therein). Figures [H [31 and H] present quadra- 
ture spectra of LMC-SCl-105 for the most prominent He I and He II lines. At phase 0.27, the 
He I lines of the secondary are significantly stronger than at phase 0.75, while, of the He II 
lines this is th e case only for He II A4200. Several mec hanisms have been proposed to cause 
this effect (see iBagnuolo et al.lll999l : iLinder et al.ll2007l . and references therein), which could, 
in principle, affect the luminosity ratio and masses derived. In the case of LMC-SCl-105, 
both lines are clearly resolved and therefore the position of the line centers, and consequently 
the radial velocities can be unambiguously measured. Studying the Struve-Sahade effect in 
eclipsing binaries could be valuable for understanding its origin, because the inclination and 
siz es of stars are addit ionally known. A thorough investigation of the effect, as undertaken 
by lLinder et al.l (120071 ). is beyond the scope of this paper. 



Following the criteria of IWalborn fc Fitzpatrickl (jl990l ). the spectral types of the primary 
and secondary are 08V and 08III-V at phase 0.75, respectively. The luminosity class of the 
secondary cannot be unambiguously determined, because of the strong emission in the He II 
A4686 line, superposed on the absorption line of the secondary. The Struve-Sahade effect 
further causes the spectral types of both stars to cha nge. At phase 0.25 the stars appear to 
have types 07V and 08.5III-V. iMassey et al.l (120001 ) have assigned a 08.5III spectral type 
to LMC-SCl-105 (or W28-22, LH 81-72) from their lower resolution spectra, likely having 
observed it near or during prima ry eclipse. According to the calibration of effective tempera- 
ture (Teff) with spectral type by lMokiem et al.l (120071 ) for the LMC, the primary should have 
TgfT = 35 — 40kK and the secondary Tefi = 30 — 35kK. A (sub)giant clas sification would mak e 
it ~ 1 — 2kK cooler than a dwarf with the same spectral class (see e.g. iMartins et al.ll2005l ). 
Figures [5] and [6] present the hydrogen and hehum lines at each quadrature and TLUSTY 
model atmospheres (see ^ for a range of Teg. The combination that best fits the data 
has Teffi = 35 ± 2.5kK and Teflf2 = 32.5 ± 2.5kK, with log{g) fixed to 3.50 and the pro- 
jected rotational velocities vsini to the synchronous values of 180 km s~^ and 140 km s~^ 
(as determined in j^ . These Tes values are consistent with the spectral type calibration of 
Mokiem et al.l (120071 ). The T^s grid step was adopted as a conservative error. Note that no 
combination of models can fit all the heliu m lines satisfactorily, in par ticular the singlet line s 
He I A4387 and A4922. IPuls et al.l J2005h point out that CMF GEN JmUier fc Mill"ei1 ll998h. 
and therefore TLUSTY (since both models are consi stent. seelBouret et al.ll2003l). predicts 



much weaker He I s inglet lines than FASTWIND (jSantolaya-Rey et al. 



1997 



Puis et al. 



20051 ). iBouret et al.l (120031 ) also state that "a simultaneous fit to all He I and He II lines is 
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almost never achieved". Computing a finer grid of TLUSTY models to better constrain tlie 
Tcff of tlie stars or using a unified model atmospfiere, sucfi as FASTWIND, was not pursued, 
because tfie determination of masses and radii is independent of the Tes- 



4. Radial Velocity Curve 

Two methods were used to measure the orbital parameters of LMC-SCl-105: two dimen- 
sio nal cross correlati o n (or TODCOR) and spectral disentangling. TODCOR was developed 
by IZucker fc Mazehl (119941 ) and can distingui sh small velocity sepa rations even more accu- 
rately than one dimensional cross correlatio n (ITonrv &: Davis Il979f) . Synthetic spectra from 
the OSTAR2002 TLUSTY non-LTE grid Jlanz fc Hubenvl boosh at half-solar metallicity 
were used for the cross correlation. The microturbulent velocity was fixed to 10 km s~^ and 
the helium abundance to the solar value, He/ H = 0.1 by number. The Teg of the models in 
the grid range from 27500-55000 K in steps of 2500 K. The surface gravity log((?), depending 
on the exact Teg value, ranges from 3.00-4.75 dex (cgs) in steps of 0.25 dex and the micro- 
turbulent velocity was fixed at 10 km s^^. The models were rotationally broadened (20-400 
km s~^ in steps of 20 km s~^) and the instrumental broadening was applied with the rotinS 
program distributed with the TLUSTY grid. 

Initially, best fit models were computed by minimizing the residuals of the sum of 2 
models shifted appropriately to the highest S/N quadrature spectrum from UT 2007 October 
24. This procedure yielded the following best fit models (Tcs, log(5'), vsini): (35000, 3.50, 
160) for the primary and (30000, 3.00, 140) for the secondary. These models were used 
as templates in TODCOR to derive initial values for the radial velocities, orbital elements 
and stellar parameters. Subsequently, vsini and \og{g) were fixed for the estimation of Teg. 
The models used for the final analysis are shown in Figures [5] and [61 The regions around 
the Balmer lines and AA4620 — 4700 were excluded from the TODCOR analysis, as the 
former are broad and the latter region contains He II A4686 emission (formed in the wind or 
interaction region), which the hydrostatic equilibrium TLUSTY models cannot reproduce. 
The resulting TODCOR velocities are given in Table [31 

An orbital fit for the systemic 7 velocity, the semi-majgr axis and the mass ratio q 
was performed with PHOEBE (version 0.31a, iPrsa fc Zwitteiil2005r). which builds on and 
enhances th e capabilities of the Wilson-Devinney program (IWilson fc Devinneylll97ll : [Wilson 
1979[ . [l990[ ). The two highest S/N spectra were assigned a weight of 2, while the lowest S/N 
spectrum was assigned a weight of 0.5; the rest were assigned weights of 1. The values 
and their formal uncertainties, found by fixing the ephemeris from the OGLE catalog, are: 
q = 0.42 ± 0.02, asini = 38.9 ± 0.5 R©, 7 = 284 ± 3 km s~^. These imply semi-amplitudes 
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of Ki = 137 ± 4 km s ^ and K2 = 326 ± 3 km s ^, and minimum masses of Mi sin^(z) = 
30.9 ± 1.0 M0 andM2sin^(z) = 13.0 ± 0.7 Mq. The rms of the fit is 12 km s'^ for the primary 
and 8 km s~^ for the secondary, which are adopted as representative values for the error in 
each radial velocity measurement. Separate 7 velocities for each star were also fit for, but 
yielded values consistent within errors with the value above, therefore were not considered 
further. 

Another accurate me thod for deriving rad ial velocities of spectrosco pic binary stars 
is spectral disentangling (jSimon fc SturmI 1 19941 ). The program KOREL (JHadraval Il995l ) 
implements the method using Fourier transforms. In order to explore the x^ parameter 
space, KOREL was run for a range of primary semi-amplitudes (100 — 160 km s~^) and mass 
ratios (g = M2/M1 : 0.2 — 0.7), using the range AA3990 — 4965, a 4 km s~^ velocity step, 
excluding the AA4620 — 4700 region. The resulting values were found to be Ki = 140^20 

in agreement with the TODCO R results. The latter hav e more realistic 



km s 



Q 






errors and were adopted in the subsequent analysis. ISouthworth fc Clausenl (120071 ) present 
the first detailed comparison of the methods developed to derive radial velocities and find 
that disentangling is the most accurate, however they applied TODCOR on a single line. 



5. Light Curve Analysis 

Besides the OGLE J— band light curve, blue and red filter light curves for LMC-SCl- 
105, roughly corresponding to V and R, are also available from the MACHO database (ID 
81.8881.21). The MACHO light curves have 331 measurements in the red and 167 in the blue 
and span 7.2 years, but are noi sier than the OGLE-I I and OGLE-HI light curve (kindly made 
available by I. Soszynski; see ISoszynski et al.ll2008l ). which has 750 points (after removing 
outliers) that span more than 11 years. The instrumental OGLE-HI light curve was offset 
to match the OGLE-II light curve in the analysis below. 

Detailed simultaneous modeling of the 3 light curves was performed with PHOEBE. 
Both the detached mode and the semi-detached mode with the secondary filling its Roche 
lobe were considered. Note, the primary star (star 1) in both binaries analyzed herein 
is defined photometrically, as the hotter star producing the deeper eclipse at phase zero. 
The parameters that were allowed to vary are: the inclination i, Tg^, the period P, and 
the surface po tentials fi^ and Q2- The time of primary eclipse Tq was fixed to the value 
determined by IWyrzykowski et al.l (120031 ) . while the value of the mass ratio and semi-major 
axis were fixed to the values determined from the o rbital fit. The light ratio was computed 
rather than fit for, following iPrsa fc Zwitterl (120051 ). Synchronous rotation, a circular orbit 
and no third light were assumed, as there was no evidence to the contrary. Tefji was fixed 



to the value determined from the spectra in ^ while the gravity darkening exponents were 
set to 1 f or stars w 'ith ra diative envelopes (i.e. T^s oc g^'"^^) and the albedo values to 0.5, 
follo wing iHilditch (J200ll ). Limb darkening coefficients for the square root law were fixed 
from IClaretl (J2000|) and the approximate refiection model was used. 



The values and errors of the parameters and the stability of the solution were explored 
using PHOEBE's scripter. The Wilson-Devinney differential corrections minimization was 
run 1000 times, each time updating the input parameters to the values determined in the 
previous iteration. The final values for the parameters were determined by calculating the 
mean and standard deviation of the resulting values from each iteration; the resolution step 
was adopted as the error for the period. The errors are similar to the conventional errors 
from Wilson-Devinney. A semi-detached configuration with the secondary filling its Roche 
lobe yielded a consistent model. Figures [3, [S] and present the phased OGLE and MACHO 
light curves and best fit model from PHOEBE. The more accurate and better sampled OGLE 
light curve exhibit s a depression before the primary eclipse, which is often seen in Algol-type 
binaries (e . g., se e iHilditch et al.l 120051 ) and is attributed to a mass-transfer stream. The 
O'Connelll (1l95ll ) effect explains the residuals near phase 0.25, since the first quadrature is 
brighter than the second by 0.01 mag. The earlier spectral type derived at this phase for the 
primary correlates with this extra fiux. All bands display deeper eclipses than predicted by 
the model. This could be due to an inac curate conversion of the differential flux light curves 
to magnitudes caused by blending (see IZebrun et al.ll200ll ). however it is also seen in the 
MACHO light curves, obtained with proflle fltting photometry. The only way to obtain a 
model with flat residuals (without changing the mass ratio) would be to additionally model 
cool spots. This was not pursued, as adding free parameters would improve the flt, but not 
signiflcantly change the masses and radii. 

Figure (TU] shows the radial velocity curve and best flt model from PHOEBE, while 
Figure [TT] presents histograms of the PHOEBE scan results and the flnal parameters, which 
are also listed in Table HI The ephemeris is: 



To(HJD) = 2450451.90113 + 4.250806(1) x E, 



'1^ 



reflning the recent period determinations of iDerekas et al.l (120071 ) and iFaccioli et al.l (120071 ) 
that are based on th e MACHO data. The period was also calculated using the analysis of 
variance algorithm of IS chwarzenb erg- Czernyl (119891 ) on the OGLE light curve and was found 
to be 4.25083 days. The uncertainties for the parameters are similar to the formal errors 
from PHOEBE. Teflf2 was left as an adjustable parameter in PHOEBE as a proxy for the 
flux ratio, however, its value and uncertainty are quoted from the spectral type calibration 
described in fj3l The light (or luminosity) ratio is L2/ Li = 0.46 ± 0.02 in J, 0.45 ± 0.04 
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in R and 0.45 ± 0.05 in V and was fixed to 0.45 in the TODCOR analysis. These formal 
errors were computed by iterating PHOEBE one r nore time with the prim ary light levels as 
additional adjustable parameters, as described by iPrsa fc Zwitterl (120051 ). The flux ratio is 
F2/F^_^0M± 0.03, in agreement with the value (0.739 ± 0.004) found for LMC-SCl-105 



by lMazeh et al.l (120061 ) in their automatic analysis of OGLE eclipsing binaries. 



The physical parameters for the system are presented in Table |5l with the final values 
for the masses and radii of Mi = 30.9 ± 1.0 Mq and Ri = 15.1 ± 0.2 Kq for the primary, 
and Ma = 13.0 ± 0.7 Mq and Ra = 11.9 ± 0.2 Rq for the secondary. The log(^) value for 
the secondary is lower than that for the primary, in support of the star being evolved with 
a luminosity class between III-V. Given the physical parameters, the synchronous rotational 
velocities of the stars are found to be 180 km s~^ and 143 km s~^. Note, the spectra were 
rotationally broadened to essentially identical values in §21 The measured masses and radii 
of the components have been determined to ~ 5% and ~ 2% accuracy, respectively. As 
discussed in §1, only 14 accurate measurements currently exist for very massive stars, with 
this work contributing towards increasing this sample to 15 (see Table [1]). 



6. Evolutionary Status 

The semi-detached configuration of LMC-SCl-105 with the less massive star filling 
its Roche lobe, along with the main sequence classification of the primary and possible 
(sub) giant classification of the secondary, point to the system being in a slow-mass-transfer 
stage of case A binary evolution. The depression in the light curve before primary eclipse 
(described in ^ also indicates the presence of a mass-transfer stream. The emission in 
the He II A4686 line and some of the Balmer lines (in particular Ha) originates near the 
secondary and can be explained by gas being transfered onto the primary. The velocities 
of these emission lines coincide with the secondary star but are shifted, indicating an origin 
near the Lagrangian LI point (see Figures [5] and [6]) . The evolutionary state of the stars can 
also be inferred by comparison with theoretical stellar evolution models. Figures W2\ and [T^ 
compare t he parameters o f the b inary with the widely used theoretical stell ar evolu t ionar y 
models by ISchaerer et al.l (119931 ) for single stars and the newer models by I Claret! (120061 ). 
respectively. The plots illustrate that these post-mass trans fer star s are o versized and over- 
luminous for their masses. According to the isochrones of IClaretl (12006! ) in a mass-radius 
diagram for single stars, the primary of LMC-SCl-105 has an age ~ 5 Myr, while the age of 
the secondary is > 10 Myr (see Figure !T3|) . which is not in agreement with them being co- 
eval. This apparent age discrepancy can be explained by mass- transfer. The ages of the stars 
according to the evolutionary models of ISchaerer et al.! (119931 ) in a luminosity-temperature 
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diagram are ~ 4 Myr for the primary and ~ 5 Myr for the secondary. 



Massev et al.l ((20001) have estimated the age and mass of LMC-SCl-105 (or LH 81-72, 
W28-22), which is located in the LH 81 association. They find a range of log(age[yr]) between 
5.55-6.77 for the highest mass unevolved stars and a turnoff mass around 70 M©. Figure [13] 
indicates an age of ~ 5 Myr or log(a g e[yr]) ~ 6.7 for the primary, in agreement with the 
value of 6.61 derived by iMassey et al.l (120001 ) for LMC-SCl-105. They also find a range of 
reddenings E{B — V) = 0.13 — 0.23 mag for LH 81 and a mean value of 0.15 mag. Late 0-stars 
have B — V ^ —0.3, therefore the measured color of LMC-SCl-105 {B — V = —0.16 mag) 
implies a reddening value ~0.14 mag. Interestingly, LH 81 contains several evol ved stars: the 



WC4 star BR 50, the WN4+0B star BR 53 and the BOI+WN star Sk -69°194 flMassev et al. 



20001 ). LMC-SCl-105 is found by these authors to have a mass of 39 Mq assuming it is single. 
This value is 26% greater than the dynamical mass of the primary or 12% smaller than the 
total mass of its components, illustrating the systematic error encountered in the study of 
massive stars, which could also affect estimations of the masses of their parent clusters. 

Given that LMC-SCl-105 has exchanged mass, a comparison with models of close binary 
evolution is needed to infer the initial masses of the components, however models at the 
metallicity of the LMC do not c urrent ly exist. According to the solar metallicity evolutionary 
models of iNelson fc EggletonI (J200ll ). LMC-SCl-105 could have had the following initial 
masses, listed as pairs of primary and secondary masses: (25.1, 22.4), (39.8, 35.5), (50.0, 44.6) 
Mq, implying an initial mass ratio of 0.90. Lower masses for the initial primary do not result 
in the measured mass ratio. These models assume conservative mass-transfer, however the 
sum of predicted inital ma sses (>4:7Mp)) i s great er than the sum of the final, measured masses 
(44.1 Mq). Furthermore, Ide Mink et al.l (120071 ) have shown that conse rvative evolution is 
not a valid assumpt ion, by comparing 50 SMC eclipsing binaries (from iHarries et al.l l2003l : 
Hilditch et al.ll2005l ) with their grid of models for a range of mass-transfer efficiencies. A 
similar grid at the metallicity of the LMC is therefore necessary to estimate the initial 
parameters of the binary. 

What will the end product of the evolut ion of LMC-SCl-105 be? The primary has a 
large fill-out ratio oi F = 0.89, as defined by iMochnacki &: Doughtyl (Il972l ). implying that 
it will quickly fill its Roche lobe when it leaves the main sequence a nd the hydrogen shel l 
bu rning phase begins. Further modeling of LMC-SCl-105, following Ide Mink et al.l (120071 ) 



or 



Petrovic et al.l (120051 ) for the mass-transfer efficiency, could indicate if the secondary will 



evolve into a Wolf-Rayet star and subsequently explode as a core-collapse supernova or a 
common envelope phase will cause the stars to merge. 
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Conclusions 



This paper presents accurate (to better than 5%) fundamental parameters of LMC-SCl- 
105, one of the brightest blue eclipsing binary stars in the LMC found by the OGLE survey. 
The aim of this work is twofold: 1) to demonstrate that targeted surveys of the brightest blue 
eclipsing binaries in nearby galaxies do indeed select very massive stars and 2) to measure 
accurate parameters for one of these rare systems. The parameters of LMC-SCl-105 were 
determined from the light curves available from the OGLE and MACHO surveys and newly 
acquired high resoluti on spectroscopy t hat t argeted quadrature phases, in part applying the 
strategy proposed by iGonzalez et al.l (120051 ) to constrain the radial velocity curve with a 
small number of spectra. The system was found to contain a very massive main sequence 
primary (30.9 ± 1.0 M©) and a possibly evolved Roche lobe-filling secondary. The spectra 
display the Struve-Sahade effect, which is present in all the He I lines, causing the spectral 
classification to change with phase, and could be related to the mass transfer occurring in the 
system. LMC-SCl-105 could further be used as a distance indicator to the LMC. However, in 
addition to accurate radii, accurate flux (i.e. effective temperatures) and extinction estimates 



are necessary for accurate distances. Eclip sing binaries have 
and independent distance s to the LMC (e.g. 



Guinan et al. 



Small Magellanic Cloud (IHarries et al.l 12003 : 



3een used to derive accurate 



1998 



and most recently to M33 (JBonanos et al.l 120061 ). 



Hilditch et al. 



Fitzpatrick et al. 2003), the 



20051), M31 (IRibaset al.ll2005h 



The accurate parameters determined herein for LMC-SCl-105 contribute valuable data 
on very massive stars, increasing the current sample of 14 very massive stars with accurate 
parameters to 15, which despite their importance remain poorly s tudied. Such dat a serve 
as an external check to resolve the "mass discrepancy" problem, as iBurkholder et al.l (119971 ) 
have shown, and to constrain stellar atmosphere, evolution and formation models. Further 
systematic studies of massive b inaries in nearby galaxies are needed to extend the sample 
of 50 SMC eclipsing binaries (IHarries et al.l l2003l : iHilditch et al.l |2005|) to higher masses 
and metallicities and populate the sparsely sampled parameter space (mass, metallicity, 
evolutionary state) with accurate measurements of their masses and radii. The method 
of targeting very massive stars in bright blue eclipsing binaries can therefore be employed 
towards this goal. 
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Fig. 1. — Mass and radius determinations of stars in eclipsing binaries, accurate to < 10% 
and complete > 30 M© from t he literature ( s ee CT and Table [1] for references). The solid 



line i s the Z=0.02 Z AMS from LSchaller et all (jl992|); the dashed line is the Z=0.008 ZAMS 
from ISchaerer et al.l (119931 ). Note the small number of measurements for stars with masses 
greater than 30 M©, all published since 2001. 
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Fig. 2. — Quadrature spectra of LMC-SCl-105, obtained with MIKE/Magellan at phases 
= 0.75 and 0.27, displaying the Struve-Sahade effect. The primary (P) and secondary (S) 
stars are labelled in the first panel. The panels display the He I A4009 and He II A4026, 
He I A4144, He II A4200, He I A4387, He I A4471 lines, respectively. Note that the lines of 
the secondary are all stronger at = 0.27; the systemic velocity is 284 km s~^; tickmarks 
correspond to 5 A in the last panel. 
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spectively. The He II A4686 emission is likely due to gas being transfered from the secondary 
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Fig. 4. — The Struve-Sahade effect in LMC-SCl-105 quadrature spectra, as in Figure El 
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respectively. Note the narrow Galactic and LMC interstellar Na I D lines at AA5890 — 95. 
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Fig. 7. — Phased OGLE /—band light curve of LMC-SCl-105. The best fit model from 
PHOEBE (solid curve) assumes a semi-detached configuration with the secondary filling its 
Roche lobe. The residuals suggest the presence of an accretion stream and hot spots (not 
modeled), arising from mass-transfer onto the primary. 
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photometry; the best fit model from PHOEBE (solid curve) is overplotted. 
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Fig. 9.— Phased MACHO b light curve of LMC-SCl-105 offset to match OGLE /-band 
photometry; the best fit model from PHOEBE (solid curve) is overplotted. 
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Fig. 10. — Radial velocity curve for LMC-SCl-105. The TODCOR measurements are shown 
as filled circles for the primary and open circles for the secondary; overplotted is the best 
fit model from PHOEBE, denoted by a solid line for the primary and a dashed line for the 
secondary. Error bars correspond to the rms of the orbital fit: 11.8 km s~^ for the primary 
and 8.4 km s~^ for the secondary. 
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Fig. 12. — Comparison of the parameters of LMC -SCl-105 with evolu tionary tracks (sohd 
hnes) and isochrones for single stars at Z=0.008 (jSchaerer et al.lll993l ). The dotted lines, 
from left to right, correspond to 1, 3, 5, 7 and 10 Myr isochrones. Both components are 
overluminous for their masses. The cooler, lower mass secondary appears older than the 
primary, indicating that the system has undergone mass-transfer. 
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Fig. 13. — Comparison of the param eters of LMC -SCl-105 with evolutionary tracks and 
isochrones for single stars at Z=0.007 (IClaretll2006l ). The dotted lines, from the bottom up, 
correspond to 1, 3, 5, 7 and 10 Myr isochrones. Single star isochrones are not compatible 
with the measured parameters for the system, which has undergone mass-transfer. 
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Table 1. Mass-Radius Data for Very Massive Stars in Eclipsing Binaries 



Eclipsing Binary Mass ctm Radius a^ 

Name- (Mq) (Mq) (Rq) (Rq) 



Reference 



F015 A 


30.0 


1.0 


16.0 


1.0 


Niemela et al C2006) 




LS1135 A 


30.0 


1.0 


12.0 


1.0 






V1007 Sco B 


30.1 


0.4 


15.3 


0.5 


Maver et al. (2008) 




VI 182 Aql A 


31.0 


0.6 


9.00 


0.18 


Maver et al. (2005) 




WR 20a B 


82.0 


5.0 


19.3 


0.3 


Rauw et al. (2004): Bonanos et al. 


(2004) 


WR 20a A 


83.0 


5.0 


19.3 


0.3 


Rauw et al. (2004): Bonanos et al. 


(2004) 


LMC-SCl-105 A 


30.9 


1.0 


15.1 


0.2 


This Work 




LMC Sk-67°105 B 


31.4 


0.7 


13.8 


0.4 


Ostrov & Laoasset (2003) 




LMC R136-42 B 


32.6 


0.1 


6.7 


0.7 


Massev et al. (2002) 




LMC HV 2241 A 


36.2 


0.7 


14.9 


0.4 


Ostrov et al. (2001) 




LMC R136-42 A 


40.3 


0.1 


7.4 


0.8 


Massev et al. (2002) 




LMC MACHO 053441.3 A 


41.2 


1.2 


9.56 


0.02 


Ostrov (2001) 




LMC Sk-67°105 A 


48.3 


0.7 


16.9 


0.4 


Ostrov & Laoasset (2003) 




LMC R136-38 A 


56.9 


0.6 


9.3 


1.0 


Massev et al. (2002) 





M33 X-7 



70.0 



6.9 



19.6 



0.9 Orggzjt al. (2007) 



*A and B denote the primary and secondary components of the EB. 
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Table 2. Log of Spectroscopic Observations 



UT 


Telescope/ 


Exp. time 


S/Nat 


Date 


Instrument 


(sec) 


4500A 


20051202 


DuPont/Echellc 


2x1200 


45 


20051204 


DuPont/Echellc 


2x1200 


60 


20051206 


DuPont/Echellc 


2x1200 


50 


20060126 


Clay/MIKE 


1x200 


30 


20060131 


DuPont/Echelle 


1x1200 


40 


20060202 


DuPont/Echellc 


1x1200 


35 


20060204 


DuPont/Echelle 


1x1200 


40 


20071024 


Clay/MIKE 


1x600 


90 


20071026 


Clay/MIKE 


1x600 


80 



Table 3. Radial Velocity Measurements 



HJD 

-2450000 



Phase 



RVi 

(km s^^) 



(O-C)i 
(km s^"'") 



RV2 

(km s^"'^) 



(O - C)2 
(km s^"'^) 



3706.76777 


0.71 


415 


3708.67764 


0.16 


185 


3710.79198 


0.65 


400 


3761.56453 


0.60 


335 


3766.55222 


0.77 


420 


3768.61519 


0.26 


140 


3770.62013 


0.73 


430 


4397.80634 


0.27 


145 


4399.84620 


0.75 


415 



1 


-25 


16 


550 


5 


20 


27 


110 


1 


-20 


-8 


600 


12 


-25 


-5 


600 


-5 


-50 





1 

4 
29 
12 
-3 

8 
-1 
-15 
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Table 4. Results From Light and Radial Velocity Curve Analysis 



Parameter 



Value 



Period, P 
Inclination, i 
Surface potential, Qi 
Light ratio in /, L2/L1 
Mass ratio, q 
Systemic velocity, 7 
Semi-major axis, a 
Semi-amplitude, Ki 
Semi-amplitude, K2 
Fill-out ratio, Fi 
Radius, ri^poie 

I"!, point 

ri^sido 

flyback 

ri'^ 



4.250806 ± 0.000001 days 

89.9 ± 0.9° 

3.073 ± 0.003 

0.45 ± 0.02 

0.42 ± 0.02 

284±3kms-i 

38.9 ± 0.5 Rq 

137±4kms-i 

326±3kms-i 

0.885 ± 0.003 

0.373 ± 0.001 

0.416 ± 0.001 

0.389 ± 0.001 

0.402 ± 0.001 

0.389 ± 0.001 



^Volume radius. 



Table 5. Physical Parameters 



Parameter 



Primary 



Secondary 



Mass (Mq) 
Radius (R©) 
log{g) (cgs) 

Teff (K) 

log(L/L0) 



30.9 ±1.0 
15.1 ±0.2 

3.57 ±0.13 
35000 ± 2500 

5.49 ±0.14 



13.0 ±0.7 
11.9 ±0.2 

3.40 ±0.20 
32500 ± 2500 

5.23 ±0.16 



